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Abstract. Strong surface magnetic fields are ubiquitously found in M-dwarfs with 
mean intensities on the order of few thousand Gauss - three orders of magnitude higher 
than the mean surface magnetic field of the Sun. These fields and their interaction with 
photospheric convection are the main source of stellar activity, which is of big interest 
to study links between parent stars and their planets. Moreover, the understanding of 
stellar magnetism, as well as the role of different dynamo-actions in particular, is im- 
possible without explaining magnetic fields in M-dwarfs. Measuring magnetic field in- 
tensities and geometries in such cool objects, however, is strongly limited to our ability 
to simulate the Zeeman effect in molecular lines. In this work, we present quantitative 
results of modelling and analysis of the magnetic fields in selected M-dwarfs in FeH 
Wing-Ford lines and strong atomic lines. Some particular FeH lines are found to be the 
excellent probes of the magnetic field. 



1. Introduction 

Magnetic fields in non-degenerate stars are found all across the Hertzsprung-Russell 
diagram, from hot high-luminous stars down to cool and ultra-cool dwarfs (see, for 
example, the review by Donati & Landstreet 2009, and references therein). The charac- 
terization of the magnetic fields in cool low-mass M-stars is of high interest because a) 
these stars contribute up to 70% by number of the total stellar population of our galaxy; 
b) these stars often show dynamo-generated activity in their atmospheres similar to that 
of the Sun (i.e. with flares seen as strong emission in X-ray, Ha, Ca II H & K lines, 
see Berdyugina (2005)), however, objects later than M3.5 (T eS « 3400K, M< 0.35M o ) 
are believed to become fully convective, therefore different dynamo mechanisms need 
to be involved to explain their fields; c) these stars are promising targets to search for 
Earth-size planets inside their habitable zones, and thus knowledge of the magnetic ac- 
tivity of a parent star is potentially important for the climate modelling of these distant 
earth's. In this paper we will present first quantitative and independent estimates of the 
surface magnetic fields in selected M-dwarfs based on synthetic spectra modeling in 
FeH lines of the Wing-Ford F 4 A - X 4 A transitions. 
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2. What is known about magnetic fields of M-dwarf up to now? 

Measurements of the magnetic fields in cool stars relies on polarimetric observations 
and analysis of Zeeman broadening of spectral lines. Utilizing the Zeeman Doppler 
Imaging (ZDI) and Least Square Deconvolution (LSD) techniques as applied to Stokes 
V spectra, Donati et al. (2008) and Morin et al. (2008) carried out a systematic char- 
acterization of the magnetic field topologies in a sample of M-dwarfs. Authors found 
a clear transition in the magnetic field topologies between partly and fully convec- 
tive stars, with later tend to host strong, poloidal, mostly axisymmetric fields contrary 
to much weaker, non-axisymmetric fields with strong toroidal component for the for- 
mer. In contrast, most recent investigation revealed a presence of M-dwarfs with strong 
toroidal non-axisymmetric fields among fully convective objects, and thus no clear tran- 
sition between these two groups of M-stars seem to exist any more (Morin et al. 2010). 
These results must be taken with certain caution, since using LSD technique and lim- 
iting only to Stokes V may result in missing an important information about the true 
magnetic field geometry and total magnetic flux, as discussed in the talk of A. Reiners 
(this meeting) and recent study by Kochukhov et al. (2010). 

The Zeeman effect is sensitive to the surface averaged magnetic field modulus 
and thus naturally gives information about the true magnetic field flux. Strong fields 
up to « 4 kG were then reported for some M-dwarfs based on the relative analysis of 
magnetically sensitive Fe I line at A = 8467 A (Johns-Krull & Valenti 1996, 2000). For 
dwarfs cooler than mid - M, molecular lines of FeH Wing-Ford F 4 A - X 4 A transitions 
around 0.99 fxm are usually utilized (Valenti et al. 2001; Reiners & Basri 2006). Some 
of these lines do show strong magnetic sensitivity, as seen, for instance, in the sunspot 
spectra (Wallace et al. 1998). 

3. Lande g-factors of FeH: roots of the problem 

Unfortunately, most lines of FeH are formed in the intermediate Hund's case, which 
theoretical description is based on certain approximations. The main problem is con- 
nected with the Born-Oppenheimer approximation, which assumes a clear separation 
between the electronic and the nuclear motion in terms of energies. This approximation 
fails for FeH, and no satisfactory description of Lande g-factors exist so far. Among 
recent improvements in understanding the Zeeman splitting in diatomic molecules one 
should mention papers by Berdyugina & Solanki (2002) and Asensio Ramos & Tru- 
jillo Bueno (2006). However, theoretical g-factors still cannot describe correctly the 
broadening of great majority of FeH lines. 

An alternative solution was then suggested by Reiners & Basri (2006, 2007), who 
estimated the (|B|/) values (/-filling factor) in a number of M-dwarfs by simple linear 
interpolation between the spectral features of two reference stars, one with zero field 
and another one with the magnetic field known from Johns-Krull & Valenti (1996), but 
the error bars of such an analysis stays high (« 1 kG). 

First semi-empirical fit to FeH lines in a sunspot spectra was presented in Afram 
et al. (2008). The authors succeeded to obtained a very good agreement with observa- 
tions for selected FeH lines and presented best-fitted polynomial g-factors of upper and 
lower levels of corresponding transitions. A little earlier, Harrison & Brown (2008) pre- 
sented measured g-factors for a number of FeH lines, but were limited to low magnetic 
/-numbers only. As an example, Fig. 1 illustrates a comparison between experimental 
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Figure 1 . Comparison between observed (Harrison & Brown 2008) and theoreti- 
cal g-factors of lower (X 4 A) and upper (F 4 A) levels of FeH with rotation quantum 
number Q = 2.5 as a function of magnetic quantum number J. Black line - pure 
Hund's case (a), red line - pure Hund's case (b), green line - intermediate between 
pure Hund's (a) and (b) cases. 



and theoretical g-factors for the lower (X 4 A) and upper (F 4 A) states of FeH transitions 
with rotational quantum number O = 2.5. It is seen that the lower state tends to be 
close to the pure Hund's case (b) and the upper state is in intermediate case, while the- 
oretical predictions always give solutions between pure (a) and (b) cases, thus failing 
to reproduce g-factors for the lower state. 

Below we attempt to provide an independent measurement from a number of mag- 
netically very sensitive molecular FeH lines, which are modelled based on the formal- 
ism described in Berdyugina & Solanki (2002). 

4. Methods 

The line list of FeH transitions and molecular constants were taken from Dulick et al. 
(2003 ) ! and corrected as described in Wende et al. (2010). The VALD (Vienna Atomic 
Line Database) was used as a source of atomic transitions (Piskunov et al. 1995; Kupka 
et al. 1999). Magnetic spectra synthesis is performed using the Synmast code (Kochukhov 
2007). The code represents an improved version of the Synthmag code described by 
Piskunov (1999). The model atmospheres are from the recent MARCS grid 2 (Gustafs- 
son et al. 2008). To compute g-factors in different Hund's cases, we implemented nu- 
merical libraries from the MZL (Molecular Zeeman Library) package originally written 
by B. Leroy (Leroy 2004), and adopted for the particular case of FeH. 

5. Results 

A sunspot spectrum is probably the only trustworthy data source to test the predicted 
Zeeman patterns of FeH lines because a) the temperature inside a spot is still hot enough 
to see strong unblended atomic and FeH lines and b) very high-resolution and S/N 
observations are available. We thus made use of an umbra spectrum from Wallace et al. 
(1998). They also derived the magnetic field intensity |B| = 2.7 kG. Our fit to the atomic 



1 http://bernath.uwaterloo.ca/FeH/ 
2 http://marcs. astro, uu.se 
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lines (mostly Fe, Ti, and Cr) in the range 9800 - 10 800 A also suggests a field intensity 
|B| = 2.7 kG (model atmosphere with T eff = 4000 K). 
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Figure 2. Comparison between observed and theoretical sunspot spectra in se- 
lected regions of FeH transitions. Blue dashed line - best-fitted g-factor from Afram 
et al. (2008), red solid line - calculations with MZL library; both with a purely ra- 
dial field of B = (2.7, 0, 0) kG. Red dash-dotted line - synthetic spectra accounted 
for the horizontal field component B = (2.5, 1,0) kG (hardly seen in the figure, co- 
incides with red solid line), brown dotted line - zero-field spectrum. Labels over 
lines indicate their central wavelengths, omegas of lower and upper states (in brack- 
ets), branch, and the /-number of the lower state (in brackets). Wavelengths are in 
vacuum. 



Using the MZL library and sunspot spectra we find that the intermediate Hund's 
case (with its present treatment in MZL) is a good approximation if (/ - lower, u - upper 
states) 

1. Q/ = 0.5 

2. Q/ or u < 2.5 and 3Y > J(J + 1) for the P and Q branches 

3. 0.1 and u = 2-5 and 5Y > J(J + 1) for the R branch. 

For the rest of transitions the assumption of Hund's case (a) for the upper level and 
Hund's case (b) for lower level provide reasonable results. 
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Figure 2 illustrates a comparison between best-fitted g-factor from Afram et al. 
(2008) and our calculation for some selected lines in the 9900 - 10000 A region. Note 
that apart from the figures shown in Afram et al. (2008), we did not make an attempt 
of correcting the theoretical spectra, i.e. no filling factors were applied. The model 
atmosphere and field intensity are the same as determined previously from the metallic 
line spectra. Still, the discrepancy is found for low omega R-branch lines like FeH 
9945 A, 9962 A, etc., which indeed show splitting closer to pure Hund's cases. Even 
though we succeeded well enough to fit the width of the observed lines with the same 
field of |B| = 2.7 kG derived previously from atomic lines, this cannot be judged to 
be more physical though until new improvements in the theoretical description of the 
intermediate case will become available. 

Having fixed g-factors for important FeH transitions, we used them to derive the 
magnetic field strength in selected M-dwarfs. As an example, here we compare model 
predictions for active M-dwarfs AD Leo and YZ CMi, for which previous attempts to 
measure its magnetic field resulted in (|B|/) = 2.9 kG and (|B|/) > 3.9 kG respectively 
(Reiners & Basri 2007). Theoretical fit to FeH lines is shown in the Fig. 3 together 
with the predictions for the non-active M4.5 star GJ 1002. In particular, magnetically 
sensitive lines such as FeH 9905 A, 9906 A, 9942 A, and 9959 A clearly point to the 
lower field modulus. That the widths of these lines are well reproduced in the sunspot 
spectra (see Fig. 2) allows us to consider them as important indicators of the mean field 
intensity. In particular, increasing |B| results in the appearance of the characteristic 
feature owing to the crossed cr-components of the two FeH lines at 9906 A. Overlaid, 
these components give rise to the absorption feature which is not seen in the observed 
spectra. Consequently, weaker fields are needed to keep these lines separated. 

Table 1 . Atmospheric parameters of investigated M-dwarfs. 



Name 


Spectral 


T e ff 


v sin i 


|B| 


m 


(PI/) 




type 


(K) 


(km/s) 


(kG) 


(kG) 










atoms 


FeH 




Sunspot 




4000 


0.0 


2.7 


2.7 


2. 7 (U 


GJ 1002 


M5.5 


3100 


2.5 










GJ 1224 


M4.5 


3200 


3.0 


«2 


1.7-2 


2. 7 (2) 


YZCmi 


M4.5 


3300 


5.0 


3-4 


3-3.5 


> 3.9< 2 > 


EV Lac 


M3.5 


3400 


1.0 


3-4 


3-3.5 


~ 3.9 <3) 


AD Leo 


M3.5 


3400 


3.0 


2-3 


2-2.5 


~ 2.9< 2 > 



|B| ra - mean surface magnetic field 
(|B|/) - results of previous investigations 

(1) - Wallace et al. (1998) 

(2) - Reiners & Basri (2007), scaled from (3) 

(3) - Johns-Krull & Valenti (1996) 

Table 1 gathers the main results of the present study including several other M- 
dwarfs which were analyzed applying the same approach. An average surface magnetic 
field resulting from the analysis of atomic and FeH lines are shown separately. For 
GJ 1224, the measurements of the magnetic field from atoms were based on Ti I lines 
located in the region 10400 - 10 800 AA, while for AD Leo, YZ CMi, and EV Lac 
results are from the Fe i 8468 A line (the same way as presented in Johns-Krull & 
Valenti (1996)). A large scatter resulting from the analysis of the Fe i 8468 A line is 
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Figure 3. Comparison between observed and theoretical spectra of selected M- 
dwarfs. Upper panel: non-active star GJ 1002, red solid line-assuming £F e = 
-4.37, blue dashed line-eFe = -4.59. Middle panel: active star AD Leo, green 
dashed - B = (1.7, 1.7,0) kG, blue dash-dotted - B = (2,0,0) kG, brawn dot- 
ted - B = (2.5, 0, 0) kG, red solid - B = (2.9, 0, 0) kG. Bottom panel: active star 
YZ CMi, green dashed - B = (0, 4, 0) kG, blue dash-dotted - B = (2.5, 2.5, 0) kG, 
red solid - B = (4, 0, 0) kG, dotted - B = (3, 0, 0) kG. 
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due to the uncertainties of fitting blue and red wings of this line, which seem to require 
different |B|. The low quality of the data (low resolution and signal-to-noise ratio) can 
also introduce uncertainties in the fitting procedure, and thus it is difficult to draw strict 
conclusions without further more accurate analysis. In addition, fitting the spectra of 
magnetic stars we assumed the magnetic field to be homogeneous everywhere on the 
stellar surface (corresponding filling factor / = 1), while previous studies assumed 
the stellar surface to consist of combinations of magnetic (f < 1) and non-magnetic 
(f = 0) regions, but with the former being represented by radial magnetic field only, 
i.e ignoring possible horizontal vector components of the magnetic field. We refer the 
interested reader to the work of Shulyak et al. (2010) for more details. 

6. Summary 

In this study we developed an approach of modelling the Zeeman splitting in FeH lines 
of Wing-Ford F 4 A - X 4 A band which then applied to measure magnetic fields in se- 
lected M-dwarfs. The main result of this study is that the magnetic field strengths 
derived from FeH lines are 15 - 30% lower than results presented in Reiners & Basri 
(2007), which are based on atomic line analysis scaled from Johns-Krull & Valenti 
(2000). The estimates of the magnetic field modulus from the Fe I 8468 A line seem to 
be systematically higher than those from FeH lines. This needs further more extensive 
investigation. 
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